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velop a neutrino oscillation mechanism induced by the presence of nondiagonal terms 
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1. Introduction 

In the development of physics, symmetry principle is a powerful tool in the construc- 
tion and interpretation of physical laws of nature. Various efforts have been dedi- 
cated to the searches of new symmetry principle beyond standard model (SM) gauge 
symmetry and ordinary Lorentz symmetry, such as SU (5) and 50(10) in grand uni- 
fied theory (GUT), or 50(32) and E{8)^E{8) in string theory. Aside from these 
gauge symmetries, ordinary Lorentz symmetry is also extended to SO{9, 1) in string 
theory or 50(10, 1) in M-theory 

On the other hand, symmetry principle is not implemented trivially in nature. 
The discovery of non-conservation of parity in 1957 ' makes people to realize that 
some sacred symmetry may be only a good approximation. And the later discovery 
of electroweak theory teaches us that symmetry could be hidden in vacuum con- 
densation, in other words, it is realized through spontaneous symmetry breaking 
mechanism. 

* Corresponding author. Email address: mabq@phy.pku.edu.cn 
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Does similar situation happen in the case of Lorentz symmetry? This is a rather 
deep question since Lorentz symmetry is a fundamental spacetime symmetry and 
has been incorporated into the two cornerstones of current physics: general relativ- 
ity and quantum field theory. The possibility to think of Lorentz symmetry breaking 
may be traced back to Diract^ through reintroducing aether into the theory of elec- 
trodynamics in the early 1950s. There are also other perspectives related to Lorentz 
symmetry violation It is first demonstrated by Kostelecky and Samuel that 

spontaneous Lorentz symmetry breaking may happen in string field theory via un- 
stable vacuum triggered by tachyon field^. After then Kostelecky and CoUaday 
incorporated Lorentz symmetry violation into the effective theory fr ame worll^J, 
which is the so called standard model extension (SME). In that wor sponta- 
neous Lorentz symmetry breaking is triggered by nonzero vacuum expectation value 
(VEV) of a tensor field in underlying theory, and these VEV of tensor fields are 
incorporated with SM fields into all possible LV operators. In addition to string mo- 
tivated LV, other approaches of quantum gravity also indicate some signatures of 
Lorentz violation. That include spin-network calculation in loop quantum gravitjEl, 
foamy structure of spacetime^, noncommutative quantum field theorjEl, and emer- 
gent gravit}EI^ etc.. However, without a complete theory of quantum gravity, all 
indications of LV above do not provide a firm and definite evidence that Lorentz 
symmetry is indeed breaking, or in other words, why it should not be an exact 
symmetry. However, we can take a positivism viewpoint that we can rely on exper- 
iments to verify or put bound on LV, as current experiments have already reached 
the sensitivity to Planck mass suppression (e.g., for dimensionless couplings, the 
sensitivity to Planck mass suppression means sensitivity to — ^ 10~^^). 

The purpose of this paper is to derive a set of modified dispersion relations 
(MDR) in the framework of SME and explore their consequence in the propagation 
properties of free particles. So we first briefly review the basic principle of SME and 
its quantum electrodynamics (QED) subset in sections 2 and 3. Then, by focusing 
on the physical relevant LV couplingJIS we derive the MDR with CPT even and 
CPT odd LV couplings respectively, together with referring their physical resem- 
blance with other distinct physical subjects in section 4. In section 5 the MDR is 
applied to neutrino propagation and ultrahigh energy cosmic ray (UHECR) prob- 
lems separately. In that section we formulate a neutrino oscillation mechanism in 
2-spinor formalism and derive a much stringent bound on the LV coupling involved 
in UHECR problem. In section 6 we give a brief summary. The convention adopted 
in this paper is 7700 = +1 for i-jf^^, and £0123 = +1 for epcrf^i,- 



2. Principle of SME 

The basic principle of standard model extension (SME) is that, SM is regarded 
as a leading order Lagrangian in the low energy effective field theory originating 
from a presumed existing fundamental theory. While the other terms are treated 
as perturbation denoting tiny departures from exact Lorentz symmetry. It is these 
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perturbation terms revealing the possible signature of physics beyond SM. And the 
whole Lagrangian could be written as 

C^CsM + 5C, (1) 

where 6C is the Lagrangian denoting tiny LV effects. Generally it has a form 

6C^^{T)^T{idfx: (2) 

where the Lorentz indices of VEV of tensor field (T) and partial differential opera- 
tors id in (2) are suppressed. These indices are matched so that they are contracted 
exactly, which indicates that SME is apparent Lorentz covariant. This is an ex- 
plicit assumption of SME (i.e. LV terms are required to be Lorentz covariant in 
their apparent Lorentz indices), and also a direct consequence of the assumption 
of spontaneous Lorentz symmetry breaking of an underlying Lorentz covariant the- 
ory, such as string field theorji^. However, this covariance property of LV operators 
should not be confused with particle Lorentz violation they indicated. According to 
the work of Kostelecky and CoUaday "'^ "'^ , observer Lorentz symmetry is nothing but 
the equivalence relation of different coordinate choice, though appropriate choice 
of coordinate system would largely simplify our calculation and in some cases even 
would be helpful in the interpretation of corresponding physical properties. While 
particle Lorentz symmetry is a real symmetry concerning the properties of identical 
particles (or localized fields) with different spin orientation and momentum through 
particle rotation or boost performed in a specified inertial frame. In ordinary theory, 
this is just the symmetry classifying different species of identical particles. While 
in LV theory, particle Lorentz transformation leaves tensor VEV {{T)) unchanged, 
thus changes the relation or interaction between SM fields with background tensor 
fields, so new phenomena may arise. In this sense particle defined as irreducible 
representation of Lorentz group is just a good approximation if particle Lorentz 
symmetry is indeed violated. 

Aside from the requirement of observer Lorentz invariance, other restrictions 
may also help us to restrict or classify LV terms. We could require the theory 
to be gauge invariant under a particular gauge group transformations, e.g., gauge 
invariance under gauge group SU[i)c ® SU{2)l®U{1)y, and that is why the 
theory is called SME. Dimensional counting may be used to classify various LV 
operators. If restricted to dimension 3 or 4 terms, this is the minimal version of 
SME or igin ally appeared i We note that dimension 5 operators are also classified 
recentlyi^l. Hermiticity and energy positivity are also necessary to make the theory 
physically meaningful. 

In addition to the above considerations, discrete symmetry transformation can 
be applied on the LV operators to classify them into CPT even and CPT odd 
classes. They form two special irreducible representations of homogeneous Lorentz 
group respectively. Using the convention of Coleman and Glashow"'^'^, the general 
irreducible representation of homogeneous Lorentz group is marked by {A, B) , where 
A,B are two angular momentum quantum numbers. So (1, 1) is identified with CPT 
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even operators and represents traceless symmetric tensor of rank 2, while (5,5) 
is identified with CPT odd operators and represents four-vector. According to the 
argument of Coleman and Glashow-*^^, the expectation value of {A, A) operator 
grows at large energy like E"^^. Thus CPT even operators dominate at high energies. 

3. QED Subset of SME 

In this section we present power-counting renormalizable QED subset satisfying all 
the requirements discussed in the above section. We can divide LV QED into pure 
photon part and fermion part. The interaction between them is included through 

covariant derivatives 'D'^=d'^ +iq A'^ . For simplicity we confine ourselves to electrons 
though the equation derived below can be applicable to more general fermions which 
are not necessarily elementary particles. The LV QED Lagrangian is 

<5^QED = <5^photon + ^-^electron, (3) 



where 



and 



in r /'cvcn i r road /r\ 

<^*-'clcctron — '^'^clcctron '^'^clcctrom I'^J 

'5>Crc"ron 3 - + ^C^ul^l^ ^ + ^d^ui^^.^ V', (6) 

5£:il^t„„ D ~a;M>r4, ~ (7) 

including those which are not directly deducible from terms compatible with elec- 
troweak structure 

<5/:cfcctro„ ^ V.^- T>ui^ + ^ig^^-^ax^ V. ^. (8) 

All the coupling coefficients c, d, e, /, g, m^, a, b, and H above are real and 
constant parameters required by the hermiticity of Lagrangian. They are related to 
VEV of tensor fields in the underlying theory. However, not all of them are physically 
observable, and some of them can be eliminated through field redefinition. This is 
the result of the fact that there is a spinor coordinate selection freedom, which 
implies that the mathematical expression of Lorentz invariant Dirac Lagrangian is 
not uniquely determinecff^. There exists an equivalent class of Dirac Lagrangian 
which are related by the fermion field redefinition of the form 

^{x)=[l + f{x,dMx), (9) 

where f{x,d) represents a general 4x4 matrix function of the coordinates and 
derivatives. For example, we can choose f{x, d) — +ia ■ x, or its finite form "^{x) = 
exp[ia • x]xix) to reproduce —a^ip^^ip from the conventional Dirac Lagrangian 

-Coirac = ^ii^l^ dfiiJ - mtlnp. (10) 
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Thus some apparent (physical irrelevant) LV couplings can be accounted for 
by field redefinition of conventional Dirac spinors, and then can be absorbed into 
redefined fields through inverse transformation. However, this field redefinition or 
field transformation works effectively only in the absence of interaction with other 
fields or interaction between spinor comp one nts due to nonlocality problem, which 
were observed by Colladay and McDonaldl^l Fortunately, since in the following we 
will focus on extended Dirac equation and its consequence, disregard photon parts 
and set the covariant derivatives Pi, into partial derivatives du, that is to consider 
only free Dirac equations, so no such obstructions will meet when performing field 
redefinition to remove some couplings. After field redefinition, the simplified ex- 
tended Dirac Lagrangian which contains only physically relevant parameters (some 
rearrangement is performed to include original Dirac Lagrangian to form a compact 
and elegant form) is written as 



/Iclectron = 9^ -^M^, (11) 



where 



= 7^ + ^"^h. + ^^"^^757. + ^5^''^^A., (12) 

M = m + 6^757^ + ^ifp.fT^". (13) 

All the coefficients in the above two equations (12) and (13) have been redefined, 
thus do not have the symmetry properties as their original ones (wi thout tilde) in 
their corresponding Lorentz indices. For details, see Colladay et al.lJ^. For simplicity 
we omit the "tilde" below and the reader should not confuse them with the original 
ones. 

Using the Euler-Lagrangian equation 

— -9^^-0 (14) 

we can obtain the extended Dirac equation below 

W7. + c^.7'^ + rfM-757'^ + ^^^'^) -(m + fe^757^ + ^i?M-tT^'')]«'(x) = 0. (15) 
In subsequent section, we will discuss this equation in detail. 



4. Derivation of Dispersion Relation of Extended Dirac Equation 

In order to get a modified dispersion relation, we proceed with the usual squaring 
procedure (which leads to Klein-Gordon equation when we apply it to the usual 
Dirac equation) to see the consequence when apply it to the extended Dirac equation 
(15). Multiplying (15) by [-iTP dp -M], we get equation 

[rT^ dpdp +i[rp, m] dp +m^] = o, (i6) 
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where 

TT^ 9p9^= {GpJg^'^p^ - d/d'^P - 2tj5<7uad''PG^^P^ - €x,ao.9^''^{G^''P^ + 75d"'')757" 
+ ^.9^"'^.9''"''[«eA^^a75 + {riaWo.u - r^^uVo^x)]} pd^„ (17) 

+ ^Hp,Hp„[zeP''p''j5 + Wn""" - n'^'n^nX. (18) 

2GPPb^-f5 + iRP^G'^P + l^d"P)['n^..lu ~ Vualf.) ~ 2id''PbP<7^,. (19) 

In equation (17) and (19) we use the definition 

[r^^P + c^''p)] = G'"'p\ (20) 

Note that, with the symmetry property of c^''p\ we have 

QP)Q(i^t^) = + c/c^"^^) + 2c'^P^'K (21) 

Obviously this squared extended Dirac equation with 2 classes of undetermined LV 
parameters ( C(^^), d^iy, H^^] b^, gxuii) is too complicated to be diagonalized in 
spinor space by continuing the same procedure. However, since their CPT properties 
are distinct, we can discuss them separately. But we should keep in mind that this 
is just for convenience. There is no fundamental reason to forbid the appearance of 
the CPT odd operators unless one imposes CPT symmetry as a custodial symmetry 
survived after the Lorentz symmetry breaking. Note that Lorentz invariance is just 
one part of sufficient conditions to the proof of CPT theorem in the local field 
theory, not a necessary one. On the other hand, CPT violation conclusi vely leads 
to Lorentz violation in local field theory, a theorem proved by Greenbergl^. So we 
could have a Lorentz violating theory with only CPT even operators, in which CPT 
odd operators are all ruled out by CPT invariance. While in the theory with CPT 
odd operators (this theory would be automatically Lorentz violating, as guaranteed 
by Greenberg's theorem ), CPT even Lorentz violating operators would be induced 
via loops, so they must be much smaller in the naive analysis with the assumption 
of tree level disappearance of CPT even Lorentz violating operators, thus can be 
neglected at the tree level calculation, which will be the case of next subsection. 

4.1. CPT Odd 

At first we write the field equation involving only CPT odd LV couplings 

• 5 + ^5am.^^^5" - ("^ + Tob ■ l)]^{x) = 0, (22) 

then by multiplying (22) on the left with [—(17 • d + ^gxpu'^^^d'^ ~ l^b • 7) — m], we 
obtain a quadratic equation 

{[Vup + |(*e^''"''75 + (^^"^7^^ - V^^^ln)9^^..9.PpW^P + b^ + 
+2j5CTp,bPd'' - ie^P^''bpgxp,^a.d''}'^ix) = 0. (23) 
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This equation can be rearranged by putting all diagonal terms (in spinor space) on 
one side, while leaving nondiagonal ones on the other side, that is 

{[^^ + \9x^..g''%^^^P + b^ + m^] + 

[^e^^"^gAM.(75 gciipdP + 47, 6^)9'' + 2-i^a p^h^ d^]}-^ {x) = 0. (24) 

Proceeding with the same squaring procedure, in principle we can get an 8th order 
differential equation without the appearance of F structure matrices. However, since 
this routine is too tedious and makes physics obscure, we do not follow this way, 
rather we concentrate on the equation (22) itself. As noted i g\^^ may arise from 
interaction among fermion constituents for a composed fermion, thus is expected to 
be suppressed further more than other LV couplings, so we simply drop it in (22) 
and get 

[i7-a- (TO + 756-7)]^'(a;) 0. (25) 
It would be easy to get a quartic order differential equation 

[(9^ + 52 _^ rri^f - 4((6 . {j,d)f + d'^)]'^ {x) ^ (26) 
from equation (25). By using the Ansatz 

"^{x) = 0(p) exp[-ip • a;], (27) 
we finally get a modified dispersion relation 

[(p2 _ 52 _ ^2-)2 _^ 4^2^2 _ . = 0. (28) 

Note that, equation (28) is noninvariant under interchange p —p, however it 
is invariant under simultaneous interchange p — s- —p and b —b. This is the 
common feature of CPT odd LV operators, which indicates a helicity dependence of 
energy levels. Taking into account of gxfj.i' just implies a further splitting of energy 
degeneracy. 

Before closing this subsection, we observe that even without LV, gravity can 
induce an equation of the same form as (25) except that, the constant vector 6^ 
is replaced by a spacetime dependent vector Ba, whose meaning will become clear 
later. This implies that gravitational field provides a practical global Lorentz sym- 
metry breaking source. The above observation will be manifested by the derivation 
of the covariant Dirac equation 

-Cspin-gravity = \/^[lV'7° " rmptp], (29) 

where the covariant derivative is Va—da —jWhcac'"^- Thus 

>Cspin-gravity = \/^[«V'7 ' ^Tp - tmplp] + [^^^j^lpj" (T'"'lpWbca] 

— drcc + jO-int , (30) 
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where £int = 4 ^ '07°'^''°V'''^fcca can be shown equal to 

^int — -Cvi + ^AI 



where W[6c]a = w^fcca and £vi is an antihermitian term, thus vanishes automaticaUy 
by the hermiticiy requirement of the theory. While £ai can be shown equal to 
^/—g'^lbld'^pB''' , where 

B" ^ {e'-'-'w^.c^a = Je'"'='e,p(aa e/ + T^^.e.'^e/). (32) 
Thus (29) could be rewritten in the form of 

spin— gravity — ^frcc H~ ^AI 

= det(e)V^[(n-a-m)+757dB'']^. (33) 

By identifying Ba with —6a, we see that the covariant Dirac equation 

[(^7 ■d-m)+ 757d^'']^ = (34) 

is of the same form as (25). 

As observed by Mohanty, Prasanna, and Lambiasel"^, the spin-gravity coupling 
can induce leptongenesis in the presence of lepton number violation interactions, 
thus may help to resolve the net baryon asymmetry problem through the so called 
electroweak sphaleron process. Since equation (34) concerning spin-gravity coupling 
and equation (25) concerning LV vector coupling are of the same form, a non- 
vanishing CPT odd LV coupling, 6^, may also be a candidate in demonstrating these 
effects and thus provide a possible solution to asymmetry problem. So the formal 
similarity of the two equations suggests that and —Ba may mimic the effects 
each other produced, thus experimental searches of the two may be complementary. 
In other words, experimental searches for curvature-spin coupling may also provide 
signals for LV bounds on and vice versa. However, there is a significant difference 
between the two. First, in LV case is treated as constant background field, while 
—Ba generated by curvature couplings is spacetime dependent though in some cases 
can be treated as semi-classical background. Second, b^ is a CPT odd LV coupling 
treated as an unaltered constant under CPT transformation, while —Ba is generated 
from gravitational sources and should transform in the same way as the ordinary 
matter field under CPT, thus makes the corresponding operator CPT invariant, 
which could be easily seen from the fact that gravitational interaction respect CPT 
symmetry. Further more, though —Ba breaks global Lorentz symmetry, it respects 
local Lorentz symmetry automatically in an appropriate free fall inertial frame, 
thus it is actually a local Lorentz invariant term. The last point is that, —Ba is 
a universal coupling reflecting a curved spacetime effect provided we insist on the 
equivalence principle to be still hold in quantum-gravity interplay region, while there 
is no good reason to regard 6^ as universal. So experimentally we could distinguish 
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the two with different physical meanings by the effects produced by nonvanishing 
Sb/j, = b'fj_ — bjj_ (where i and j refer to different flavors). 

4.2. CPT Even 

Next we turn to CPT even LV corrections to conventional Dirac equation and derive 

the corresponding dispersion relations. 

For completeness, we rewrite the CPT even Dirac Lagrangian below 

-CSron = ^i^iv^^u + c^u ' d^uln)^ 9^ V " ?^(m + H^. cj^i', (35) 

where we replace H/^^ by H^^,,, i.e. 

H^^u= H^^ X M. (36) 

The meaning of this replacement will be clear later. Then we write down the corre- 
sponding field equation by setting gx^i, and b^ equal to zero in (15), that is 

[ih. + c^.r + d^ulni^) -(m + hl^,u^")\^{x) = 0. (37) 

Multiplying (37) on the left with [^(7^ + c^^i'' + d^^TsT'') -^iJ^^cr^'^ + m], we 
get the following equation 

{-[Gp^G", - dp^d", - 2mp„d%Yl^G% + d%Wd^ + %eo.paAl^G\ + dW^^^W'^d'' 

+ -^i^ocpa^l^H'^'H'^^ + 27?^] _ m2}*(:c) = 0. (38) 

In deriving this equation, we use the anti-commutation relations 

{iv, (J pa) = hu, ^bpla - Icrlp]} = -^C^paa^b^", (39) 
Wpa, (^ixu) = 'iiepap.ulh + "^{VpfiVai^ " 'npuflan), (40) 

which can be proven by direct calculation. Note that, for aesthetical consideration, 

we have retained the term -2iapadP^d%d^'d'' in -2icrp^d''^[75(c + r])"^ + (^^]^^'^" , 
which vanishes automatically for antisymmetric properties of Upa- In (38) we use 
G^v defined in (20) instead of c^i^, and this definition is triggered by the observation 
that the Minkowvsky metric 77^,^ is always followed by c^,^. 

Eq. (38) shows that, it is hard to be diagonalizcd in 4-spinor space by the squar- 
ing procedure we used before. This is due to the entanglement between nondiagonal 
terms involving d^i^ and So we can derive the modified dispersion relation by 
assuming = or if^y = respectively. 

When Hfjti, = in (38), it leads to equation 

{[{rif,u + Cp^cP, + 2c^, - dp^dP,)d''d'' + m2][(?7„^ + c^ac"^ + 2c„0 - d„c.d%)d'^d^ + m'] + 
^[dpad'^ic + vTAc + VU - d'^{c + r,)p0d\{c + r,)^,]d''d''d''d^}^{x) = 0. (41) 
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This remains to be a complicated equation. By using Ansatz (27) we can get a 
quartic-order dispersion relation 

+A[dp^dP^G%G„[i - dP^Gppd\G^,]p^' p" / = 0, (42) 

and this is still a complicated equation. We can analyze the role of dun and c^i/ 
separately in (42) by setting the opposite term equal to zero respectively. 
Setting dfjtv = corresponds to 

[{ri^.u + Cp^cf, + 2c^,)p^ - = (G^^G^.p'^ - m^) 

= {G^np'' p" -m^)=Q, (43) 

where at the last step we define G^j/ = Gp^G''^. This definition makes (43) looking 
like a formally rclativistic dispersion relation, except with Minkowvsky metric 77^1/ 
replaced by G^„. From (43) and the definition (20) we see that c^i, behaves like a 
fluctuation of metric. 

While setting c^i, = in (42) corresponds to 

[(p2 + ^2 ^ dp^dP,p^^ p'^f - 4(m2p2 + (dp^pP pP f)] = 0, (44) 

which is a quartic order equation, but could be solved formally by using triangle 
parametrization with 

= 4m V, 

= 4idppp" p>')\ 

(p2+m2 + rf^^rf^^p''p-)2. (45) 

Thus + = Z^, which is just the identity of (44), and set 

y = Zsin[6l], X = Zco%\d\ (46) 

with B ^ 1 since dp,j constrained by experiment must be small. 
When dfj,v = in (38), it leads to equation 

[e„pafse^Cr>5H"fH^'^{^H-^H^' - rj'"' G^G^d-" 8^)]}^ (x) = 0, (47) 

where = Hi^iH^r,- 

When replacing id^ — > in (47), we can solve the corresponding quartic order 
equation in the quadratic form 

+ (cppc^, + 2c^,)p'^ p" + ^H^ - = 0, (48) 
ec^p^pe^^r^sH'^PH'^^i^H^fH^' + r,^' G'^^G^p" p^) = 0, (49) 



September 4, 2009 11:32 WSPC/INSTRUCTION FILE Iv-final 



11 

where (49) should be interpreted as a constraint equation on H. It is possible to 
separate (48) from (47) only in the case of the assumption that (49) is semi-positive 
defined. This assumption is satisfied definitely when the vector defined below 

= ]^e^p„pH^PG%p'' (50) 

is a timelike or lightlike vector, which can be seen from the expression (49). In 
equation (49) we see H appears with a totally antisymmetric tensor e^fQriS, its ap- 
pearance reminisces us the anomaly expre ssio n in the presence of gauge field, that is 
A{x) — —j^efiypaFj^'^{x)Fp''{x)tr[tatfit^P^. So we guess the physical effects pro- 
duced by if^'^ defined in (36) (rather than H^'^) resemble that of electromagnetic 
field. This is indeed confirmed by the operator — jiyV' Hiiv u'^'^tjj in the Lagrangian 
(35), which is nothing but the "Pauli term"Il^ appeared in the Lagrangian form. This 
term could give an additional contribution to fermion magnetic moment (anomalous 
magnetic moment), thus could be constrained by muon "g-2" experiment. 



5. Applications 

Since we have already derived a set of modified dispersion relations induced by 
various LV couplings, we can see what novel physical consequences these relations 
can lead to. It is well known that = ~p'^ + is a fundamental equation in 
conventional physics, so modification of this equation is expected to have a wide 
impact on high energy physics at which LV effects are expected to be less suppressed 
than at low energies. Indeed, introducing even minuscule LV would lead to processes 
conventionally forbidden at high energies, or accumulating unexpected observable 
effects when particles propagate through cosmological distance, or even lead to some 
processes allowed at intermediate energy range while forbid den at higher and lo wer 
energieJlQl. For example, radiative muon decay /i — > e + T^^, neutron stabilitjMl 
at ultrahigh energies, and va cuu m dispersion and birefringencJ^, vacuum photon 
sphtting-'^^, and photon decajff^, etc. 

However, for the dispersion relations we derived, there are two points to be clar- 
ified. Firstly, there is a preferred frame in which each relation has a most simplified 
form. So when we use this form of relation, we implicitly presume that a preferred 
inertial frame has been chosen. Secondly, as we previously commented, the LV cou- 
plings involved are redefined couplings, while this redefinition only works properly 
in the absence of interactions. So, strictly speaking, these relations are only ap- 
plicable to propagating problems where particles involved could be considered as 
free fermions, though no obstruction would meet in deriving a dispersion relation 
from the interaction equation, where partial derivatives are replaced by covariant 
derivatives and radiative corrections are include to calculate a complete propagator. 
Then, since CPT-odd equation resembles the covariant Dirac equation with spin- 
gravity coupling involved, which had been extensively discussed elsewher4^, and 
CPT-odd operator grows with energy increase much slower than CPT-even one. 
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which has been discussed in the end of section II, we focus our attention on CPT- 
even couphngs and discuss the imphcation of the corresponding MDR in neutrino 
and ultrahigh energy cosmic ray problems. 

5.1. Neutrino Oscillation 

Neutrino oscillation might be the only definite signal indicating physics beyond 
SM, and has been extensively discussed in the literaturJlSI, While most of them 
focus on neutrino with mass nondegenerate scenario and use Dirac equation as a 
starting point, which makes the assumption of neutrino classification unclear. In our 
derivation, we assume neutrinos to be Majorana spinors from the beginning, thus 
it is suited to be described in a 2-spinor formalism. As this assumption indicates, 
we should first reduce the Dirac equation into a 2-spinor form. Beginning with (37) 
by ignoring H^^ for simplicity, we rewrite the equation as 

[«(7. + c,,.7^ + df^.Toin -m]^{x) = 0. (51) 

For convenience, we redefine these LV couplings in a manifest V-A form resembling 
the V-A theory, which is a low energy effective field theory of electroweak theory. 
The definition is 

gtu = (c - d)^,, = (c + d)^,. (52) 

With this definition, equation (51) could be written as 

[^{l•^+9L^-^r + 9'',^-^r) 5^ -"^]*(^) = 0. (53) 



Using the projection operator and the definition of Majorana spinor ^! = ^'^ = 
T 

, we can derive from (53) the corresponding equation 

+ Cf,^(T^') (T (j)- ima^cf)* ^ (54) 

satisfied by Majorana 2-spinors, where = (— 1, cr) and (j> is the reduced wave 
function. For details, see Appendix. Rewrite (54) in the form of Schrodinger equation 

d 1 ^ ^ 
i—(f> = - a - V (t>- ici.^a^' d" 0, (55) 
o t I 

where we already assumed neutrino to be massless fermion, as our derivation of 
neutrino oscillation will not be based on mass nondegenerate scenario. Note that 
is a simple notation of a column of 2-spinors, and in our case, we consider only 



two flavors as an illustration, so (j) = [ \ ^ ] , and the corresponding Hamiltonian 

H — ^ a ■ —ic^iyd^^d'^ should be regarded as a 2 x 2 matrix operator in flavor 
space. Using the Ansatz (j){x) = 4>{p) exp[— zp - a;], we can write the Hamiltonian in 
momentum space as 

V -{c^.u)2l<T^' P" ^ -P -M22<y^ J ' 
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We can diagonalize this Hamiltonian by a rotation matrix R. Choose a spe- 
cific reference frame in which rotation invariance still holds, then we can assume 
{cij)a[3 = kafjSij (where i and j run over 1 to 3, a and (3 run over 1 to 2), and that 
all the other terms are zero. Since Lorentz violation is stringently restricted to be 
tiny, kai3 «C 1, we can simply drop the diagonal terms of {c^v)ai3 in (56), and assume 
fci2 = ^21 — k. Then we can get the corresponding eigenvalues of Hamiltonian (56) 
as 

Ai = l + fc, A2 = l-fc, (57) 
and the corresponding rotation matrix 



8m[9] cos[e] J x/2 Vl 1 

which just corresponds to set the rotation angle in (58) as 6* = ^. Then we can get 
the relation between energy eigenvector R(j) with 



cos[6'] sin[6i]\ f\u)i 



Vt) J \ — sin[0] cos[6'] J \ \iy)2 



(59) 



For a muon-type neutrino emitted, the neutrino state after evolving through a time 
interval t is determined by 

\v{t)) = {siii[9]exp[~i5E ■ t]\i^)2 + cos[0]\iy)i) exp[-iEit] 

= sin[26'](exp[-i(5£; • t] - l)\i^r) 

+ (1 + sin[6l]2(exp[-i(5£: • t] - l))|i/^)] eyip[-iEit], (60) 

so the flavor transition probability is given by 

- sin^[20]sin^[^&^], (61) 

where at the last step we used (57) and &E = E2 — E1 = —2k\~])\. Including the mass 

terms in (55) and (56) just complicates our formula without any principled difficulty. 

Note that the presence of nondiagonal tensor couplings ((c^j/)q/3 ^ for a ^ (3) 

in flavor space is essential for this Lorentz violation induced neutrino oscillation 

scenario. Though this is a rather simple model to illustrate neutrino oscillation 

caused by tiny Lorentz violation, we can still gain some insight by comparing it 

with experiments. The MINOS experiments reported their oscillation fit results with 

sin^ [26*23] > 0.84 with 90% confidence level and Am^g = 2.38toTe x lO^^eV^ with 

68% confidence levePSland they analyzed the data with the same two fiavor i/^ — > 

oscillation assumption. We find that the oscillation angle 6* = |- is consistent with 

MINOS results, even very close to it. From the squared mass difference (- lO^^eV^) 

Fool 

and the robust bound from cosmology on the sum of neutrino mass^ (0.5 — 1.0 eV), 
we can estimate the neutrino mass to be around 0.1 eV order. Then Lorentz violation 
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coupling would contribute to an effective mass term as can be seen from (54) when 
neutrino is significantly energetic for k ■ ^ Wneutrino- In the MINOS neutrino 
experiment, the peak in neutrino ener gy spectrum is around 3 GeV (low-energy 
beam) to 7.8 GeV (high-energy beamj^O', so we can give a rather rough bound 
on the size of Lorentz violation coupling k ~ '"""-gp" ^ f^IcV 10~^°, if Lorentz 
violation contribution to neutrino oscillation in MINOS experiment is comparable to 
the non-degenerate mass contribution. In principle, this bound could be restricted to 
more stringent accuracy of order k ^ 10^^^ ^ '^'^P by dimensional analysis. Since 
this estimate of the size of Lorentz violation coupling is just the inverse of 7(= ^) 
factor of the high energy neutrino and depends on the assumption of comparably 
contribution of Lorentz violation, to obtain more accurate estimate of the order of 
this LV coefficient we need to take into account the mass effect (i.e., (54) is used) and 
give a more reasonable weight on Lorentz violation contribution by using experiment 
data (e.g., Am^g) directly or even by matching the whole energy spectrum. We 
noticed that some more comprehensive worliSQi have already been done in the three 
flavor case which involves the whole renormalizable LV operators (c, d, a, 6, e, /, g, 
and -ff ) , though this was done under some reasonable perturbative expansion (since 
Lorentz violation correction would be tiny) in order to get the required effec tive 
Hamiltonian which controls neutrino propagation effects. In their first papeJ^^l^ 
a general framework in the LV induced neutrino oscillation was given and some 
definite signals in experimental searches for Lorentz violation in neutrino sector 
were classified and examined. Their subsequent works focused on particular models 
where the number of nonzero LV parameters were reduced significantlj^^. In this 
sense, our work is just a illustration or toy-model, however, its simplicity makes the 
oscillation mechanism induced by tiny Lorentz violation more obviously and the 
assumption of the neutrino property (Majorana neutrino) more apparently. Further 
it could be generalized to more practical model (3-flavors) directly by taking into 
account mass terms since the pure Lorentz violation (i.e., massless neutrino case) 
model may not be a practical solution to globally fit al l ne utrino oscillation data 
from solar, reactor and atmosphere neutrino experiments!^. However, whether the 
generalized form could accommodate with the experiment data or not still lacks 
checking. 

Some remarks should be said about the transition probability (61) which is pro- 
portional to neutrino energy, instead of inverse proportional to it as in the case 
of mass nondegenerate scenario. This property is the common feature of all non- 
standard neutrino oscillation scenarios a nd reminisces us the neutrino oscillation 
induced by equivalent principle violatiorP^. The formula of which is exactly the 
same as (61), except replacing 2fc with where /iqo = — 2(/> = is the 

00-component of metric fluctuation, and a is the post Newtonian parameter (in 
general relativity, a = 1 and is universal). This similarity is not an accident, since 
Lorentz violation is assumed as an remnant of quantum gravity, and in quantum 
region there is some indication that equivalence principle is violated. We guess that 
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the equivalence principle violation may indicate Lorentz violation at least locally, 
as the equivalence principle ensures the existence of local inertial frames, which is 
the premise of local Lorentz transformation. Actually, Lorentz violat ion must be 
followed by equivalence principle violation, an issue recently clarified irPSl. Further- 
more, we can see from (52) that nonvanishing df^, gives rise to different couplings to 
left and right handed Dirac fermions, so it may induce observable effects in energy 
splitting between different helicities. 



5.2. Ultrahigh Energy Cosmic Ray 

Ultrahigh energy cosmic ray (UHECR) provides a natural source of high energy 
particles, with energies up to 10^^ eV, much higher than that of energetic particles 
generated by man-made accelerator. But the energy of UHECR reached earth can 
not be much higher than that, as it is predicted to be terminated at around 5 x 
10^^ eV for the energy lose in the collision of UHECR particle with CMB photons 
by Greisen, Zatsepin and Kuz'mirPS, which is known as the GZK cutoff. Similar 
situation happens in the collision of TeV 7 ray with infrared photons. However, 
this prediction has not yet been confirmed by experiments. AGASA, Fly's Eye both 
claimed that they observed events with energies nearer or above this cutofP^, while 
HiRei26l 

and Pierre AugeJ^H experiments recently claimed the observation of the 
cutoff. This unsettled problem has stimulated many attempts to resolve it, including 
active galactic nuclei (which is favored by Pierre Auger experiment), primary flux 
of magnetic monopoles, "Z-boson bursts" produced by collision of UHE neutrino 
with relic neutrino nearby, pseudo-complex extension of standard field theorjl22l^ 
etc.. Of course, LV also provides a possible candidate to extend or entirely rule 
out this cutoffl^. In this paper, we follow the general analysis of Coleman and 
Glasho^^ffS and show that the LV coupling in modified dispersion relation (43) 
could be constrained either from the absence or the presence of GZK cutoff. 

We take the common assumption that the identity of UHECR are protons, and 
analyze the pion-nucleon resonance formation reaction P + 7(CMB) — > A 
(1232), which is the dominant contribution to GZK cutoff. This reaction is possible 
if and only if Eq > -Emin(Po), where Eq is the total energy of initial particles 
and ii^min( P 0) denotes the minimum total energy of the final products, whose total 
momentum is equal to initial total momentum Pq, which is implicitly assumed 
from energy momentum conservation. So the reaction is kincmatically allowed by 
the condition 

uj + Ep> Ea, (62) 

where uj is the energy of CMB photon, Ep and i?A arc the energies of proton and A 
resonance, with the subscripts denoting proton and A resonance respectively. We 
rewrite dispersion relation (43) in a relativistic form 



1:12 ~E)2 2 I 2 4 



(63) 
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where Ca is the maximal attainable velocity for the ath particle defined irffQ^. In the 
case we considered, it just requires the definition below 

1 



1 + cgo 



(64) 



where Cqq denotes the 00 component of c^,y for the ath fermion, and other compo- 
nents of which are assumed to be zero. We rewrite (62) in terms of Ep and u, 



^{\Pp\-^fc\+mlc%. (65) 



For threshold reaction, the initial proton momentum P p is coUinear with that of 
A resonance while anticollinear with that of CMB photon, which has already been 
used in (65), i.e., the substitution of Pa = P p — lo. Squaring (65), we have 

El{l + f^)(l - £^) + 2Lo{Ep + \Pp\cl) + u\l cl) 

Cp Cp 

+ {ml - m\—)clc\ > 0. (66) 

Cp 

For ultrahigh energy proton: Ep ^ \Pp\- As LV coupling should be much smaller 
than 1, we take the approximation that ca ~ 1 and 1 + — ~ 2. Substituting this 
approximation into (66) leads to a quadratic order inequality of Ep 

2El{l-—)+'iLuEp + K>Q, (67) 



where 



(.2 

K={ml-m\^)clei. (68) 



Thus the threshold energy that the reaction kinematically allows is the small positive 
root of (67) when equality is hold. With the assumption 

1 - — > 0, (69) 

Cp 

and quadratic order equation (67) (with equality hold), we obtain two roots with 
opposite signs. One is negative but with larger absolute value, the other is positive 
which gives the threshold energy, i.e., pion-nucleon resonance formation reaction is 
kinematically allowed only for energy above this positive value. On the other hand, 
with the assumption 

1 - — < 0, (70) 

Cp 

we obtain two positive values. The small one is the threshold energy, while the larger 
one is the terminating energy of this reaction, which means that the formation 
reaction is kinematically allowed in an intermediate energy band. This is a striking 
feature of Lorentz violation corrections to the famili ar particle reaction process 
previously referred and was systematically discussed i so GZK problem in this 
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case provides one concrete illustration of the analysis irff^. However, both cases give 
the same threshold formula, that is 



1-1/25(1 



P 1- — 

4a; SZuj-^ Cp 

where the first term at the last step is the conventional threshold energy in the 
absence of LV, an d th e other terms are tiny LV corrections. Substituting -Ethic = 
— ^ = 5 X 10^^ eVf^ll and the experimentally observed threshold energy i?p = 5.6 x 
into (71), we can deduce the bound 



CA _ 2w{ Ep - £^tlirc) 



1-:^ = - """^ (72) 

'^thrc 



Since A resonance is not a spin-1/2 fermion, we simply assume ca — 1, and sub- 
stituting (64) into (72) yields the bound on LV coupling Cqq ^ 10^^^, which is 
more stringent than that derived inl*^ by two orders of magnitude. This is consis- 
tent with our expectation since we adopt the data confirming GZK cutoff. Note 
that this bound is of importance only in the sense of order of magnitude, since it 
cannot be fixed firmly from the location of GZK cutoff alone, many other effects 
could compensate for some amount of Lorentz violation, for example, the uncer- 
tainty of source distribution. However, the bound we obtained is strict enough and 
has already reached the Planck mass suppression sensitivitylO"^'^, which indicates 
that Lorentz violation in dimension 4 operators (for proton-LV tensor coupling) is 
indeed too minuscule to be detected. 



6. Summary 

Searching for Lorentz violation (LV) experimentally or theoretically has received 
much attention in recent years. As QED has been tested to a marvelous accuracy, it 
is expected to be an ideal research area to probe the minuscule trace of LV both in 
theory and in experiment. In this paper we studied various modified dispersion re- 
lations (MDR) derived from extended QED by the assumption that a particular set 
of LV couplings is nonzero. This is a reasonable assumption since if a fundamental 
theory really violates Lorentz symmetry, the corresponding tensor coupling in low 
energy effective theory should contain less parameters than what we presented here. 
On focusing these LV couplings in MDR or extended Dirac equation, we observe 
the similarity of some LV couplings with the spin-gravity couplings or metric cou- 
plings in covariant Dirac equation. This resemblance may indicate a deep physical 
relevance of these LV couplings with that in the quantum-gravity interplay region, 
since Lorentz violation is assumed as a remnant signal of quantum gravity. 

In addition to that, this similarity has also been observed in the application to 
the neutrino sector. We also derived an oscillation formalism by explicitly assuming 
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neutrino as Majorana spinor. We found that the nondiagonal terms of LV couphngs 
in flavor space involved in neutrino sector could explain neutrino oscillation even in 
the massless c ase. T hough this possibility has been comprehensively discussed by 
several 

authorJSnH] 

and even have already been partly tested in some experiment 
such as it is the first time, as far as we known, to derive the oscillation 

by explicitly assuming neutrino as Majorana spinor in 2-spinor formalism, thus our 
model could be viewed as a toy-model simply demonstrating neutrino oscillation 
induced by Lorentz violation. We also made a rough estimate on the size of LV 
couplings involved as fc 10~^°. We note that the LV couplings for different species 
of fermions involved in specific problems are defined independently since we can- 
not calculate them from an underlying concrete model displaying this LV effective 
Lagrangian as the low energy limit after spontaneous Lorentz symmetry breaking. 
In principle this oscillation formalism can be generalized to the 3 flavor case by 
taking into account the mass effect. Then by comparing it with the neutrino energy 
spectrum obtained experimentally rather than only with the mixing angle and mass 
square difference, we expect that more accurate bounds on LV couplings could be 
obtained in neutrino sector. 

By application of MDR to GZK problem, we derived a much stringent bound on 
the order of magnitude of LV parameters of protons (cnn^ 10~^^) from the recent 
observation in HiReJ^ and Pierre Auger experiments^ and we note that more 
bounds could be obtained on neutrons by taking account of some exotic process 
(such as proton vacuum Cerenkov radiation) in the analysis of the propagation of 
ultrahigh energy cosmic rays when Lorentz violation is presenlPSl. There are also 
many stringent bounds on the magnitude of various LV parameteri^S] up to date. 
We observe that most of these bounds are actually the bounds on the difference of 
LV parameters to different species of par ticles, as in the case of nonuniversal gravity 
coupling ind uce d neutrino oscillation - or the maximal attainable velocity analysis 
developed irffSl. They just indicate (from the opposite side) that the difference of 
Lorentz violation tensor couplings to different species of particles is rather small, in 
other words, they strongly suggest that the tensor fleld triggering Lorentz violation 
in the underlying theory couples to the standard model field universally, at least 
for dimension 3/4 operators. So whether Lorentz symmetry is just a perfectly good 
approximate symmetry is still an open question. 
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Appendix A. 

To get an equation describing Majorana 2-spinor in (54), we begin with (53) de- 
scribing Dirac 4-spinor. Actually, we could use 2-spinor or 4-spinor formalism to 
describe either Dirac or Majorana fermions. However, for a Majorana fcrmion the 
independent degrees of freedom are 2, thus not all of its 4 components in 4-spinor 
formalism are independent, while for a Dirac fermion, it needs at least 2 different 
2-spinors. So for a Majorana fermion, it is adequate to be described by 2-spinor for- 
malism, while 4-spinor formalism is suitable to describe a Dirac fcrmion. In order 
to get the suited equation describing Majorana 2-spinor from (53), wc need to use 
the projection operator and the Majorana spinor definition = ^f"^ = C"^ to 
project equation (53) from 4x4 matrix space to the irreducible 2x2 subspace. 
First, we give the projected wave function and F matrices below 

*^ = i^^, ^«=1±^^; (A.l) 

L _ 1-75 R_ 1 + 75 „^ 

7 = 2 ^ = 2 ^' ^ ' 

Using these definitions to rewrite (53) in the form 

[i(7? + af^.i^'n + i(7^ + gjl.j'-n - m(*^ + *^)] = 0, (A.3) 

and multiplying (A.3) from the left with 75, we have 

W7.^ + 9j!.l''n d" - ^(7^ + .9^.7''^) d^' - m(vI/« - v&^)] = 0. (A.4) 
From (A.3) and (A.4) we can get two independent equations 

W7^ + g^.l^'n - m*""] = 0, 

M + 5^.7'"') - m*^] = (A.5) 

for left handed and right handed fermions respectively, with the mass term mixing 

each other. The two above equations are not independent from each other for a 

Majorana fermion. We can take complex conjugate operation on (53) and charge 

J" 

conjugate operation = on wavefunction to get a charge conjugate equation 
of the same fermion field 

[i(7. + 9^.^^r + gfl^^^in -mW{x) = 0, (A.6) 



which is just (53) with and interchanged. Since the charge conjugate field for 
a Majorana field is just the original field multiplied with a phase factor (see (A. 7)), 
the field equation satisfied for the charge conjugate field should be the same, thus 
impose the condition g^,^ = g^^, i.e., = g^^ = c^^. Then we can use the same 
procedure above in getting (A.5) to get a set of corresponding equations for (A.6), 
which is nothing but the same equations of (A.5). Note that for Majorana field, 
equation (A.5) contains actually just one independent equation, the lower one is 
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the equivalent form of the upper one. Using the condition satisfied by Majorana 
fermion below 

^{x) = 7j^'''{x), (A.7) 

where 77 is a phase factor (for simplicity, we take it equal to 1), we can get the 
relation 

\i,R{x) = ^^^^^%x) = -7°C*^*(x). (A.8) 
Then substituting (A.8) into (A. 5), we get the equation below 

i{v,^ + c^ui") (l> + m<i>* = 0, (A.9) 

where 

r^^^ = C-S°7^'' (A.IO) 

(l){x) = ^^{x). (A.ll) 

Since (A.9) is expressed in 4 x 4 matrix space, thus could be reducible. We 
can find an irreducible representation of ry'' in 2 x 2 matrix space, that is setting 
T]^ = ia'^a^, where = {—l,~a). So in 2-spinor irreducible space, (A.9) becomes 

i(o-^ + c^uT^) d" 4)- ima'^(f>* = 0, (A.12) 

which is just (54). 

Wc can derive this equation satisfied by Majorana fermion from a more manifest 
way in displaying its Majorana feature, i.e., rewriting the modified Dirac Lagrangian 
in 2-component formalism and using the neutrality condition to drop the coupling 
terms between the two fermions. We first give the 2-component formalism of mod- 
ified Dirac Lagrangian, which is obtained by expressing 

in the Lagrangian (35), where a'^ = {—1,'a), = {—l,—'a). The resulting La- 
grangian (only contains Lorentz violating couplings) is 

/^Dirac = -ir]{c + d)^t,a^d"rj - ix{c - d)^i,a^d''x 

+ iH^4T,a''''x + W"'rj). (A. 14) 

Note a'^'^ appears here is defined to be a^" = \{a^a^ — a^a^), where a^" = 
ji'^^cr" — a'^a^), rather then that appeared in the text, which could be relabeled 

by T,!^" = |[7'*,7''] = 2i ( ^ j. For any Dirac spinor tpjj, we can construct 

two corresponding Majorana spinors as 

1 -i 
^Mi = -^{i^D + ^D), V'M2 = ■^(V'r' - V'd), (A.15) 
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T 

where iplj = CiJjD is just the charge conjugate fermion field. Using this equation 
we can decompose the two Weyl-spinors as two decoupled Majorana spinors in the 
absence of additional inter-couplings (e.g., LV-couplings) by the deduced equation 
below 

X = ■^(V'l +«V'2), ??= -^(V-i -«^2), (A.16) 
where l-Mi =( ^ 1 , i = l,2. 



Then the corresponding Lagrangian (not the full one (35)) is 

2 

+ H^A-4'2(r'"'i'i + i^a^'^'lh), (A.17) 

where we have used the fact ipa'^'^ip = 0. So we can explicitly see the couplings 
between two Weyl-spinors arising from d^i, and -ff^j/. In the Majorana theory, 
we can simply drop them. Thus the corresponding full Majorana-Lagrangian in 
2-component theory is 

-CLjorana = KV^^-^ + C^u)d'' i)a^ii + ^(V'V + V^), (A.18) 

and the corresponding 4-component form is 

^Majorana = ^*(7^ + C^'l.) 9,. * " ^ (A.19) 

From the Majorana Lagrangian (A.18), we can deduce an equation 

i{ri^i. + c^,^)a^'d"^ - = 0. (A.20) 

Prom the definition = X^m, where ^ = > ^"^^ the convention A = 1 being 

chosen (which could be confirmed from equation (A. 15)), we can deduce 

^ = ia^ii*. (A.21) 
Substituting this equation to (A.20), we again obtain equation (A. 12). 
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